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Abstract. The current status of neutrino cosmology is reviewed, from the question 
of neutrino decoupling and the presence of sterile neutrinos to the effects of neutrinos 
on the cosmic microwave background and large scale structure. Particular emphasis is 
put on cosmological neutrino mass measurements. 
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1. Introduction 

Next to photons neutrinos are the most abundant particles in the universe. This means 
they have a profound impact on many different aspects of cosmology, from the question 
of leptogenesis in the very early universe, over big bang nucleosynthesis, to late time 
structure formation. In the present review I focus mainly on late-time aspects of neutrino 
cosmology, and particularly on issues relevant to cosmological bounds on the neutrino 
mass. 

The absolute value of neutrino masses are very difficult to measure experimentally. 
On the other hand, mass differences between neutrino mass eigenstates, (mi, 777.2,777.3), 
can be measured in neutrino oscillation experiments. 

The combination of all currently available data suggests two important mass 
differences in the neutrino mass hierarchy. The solar mass difference of b~m\ 2 ~ 7x 1CT 5 
eV 2 and the atmospheric mass difference 5m| 3 ~ 2.6 x 10~ 3 eV 2 [1-3]. 

In the simplest case where neutrino masses are hierarchical these results suggest 
that mi ~ 0, m 2 ~ Sm so i aT , and m 3 ~ 5777 atmosphcric . If the hierarchy is inverted [4-9] one 
instead finds m 3 ~ 0, m 2 ~ 5777 atm os P hcric, and m x ~ <5m atm0 sphcric- However, it is also 
possible that neutrino masses are degenerate [10-20], m\ ~ 7772 ~ 7773 ^> 5777 atm0 sphcric, 
in which case oscillation experiments are not useful for determining the absolute mass 
scale. 

Experiments which rely on kinematical effects of the neutrino mass offer the 
strongest probe of this overall mass scale. Tritium decay measurements have been 
able to put an upper limit on the electron neutrino mass of 2.3 eV (95% conf.) [21]. 
However, cosmology at present yields an much stronger limit which is also based on the 
kinematics of neutrino mass. 

Very interestingly there is also a claim of direct detection of neutrinoless double 
beta decay in the Heidelberg-Moscow experiment [22,23], corresponding to an effective 
neutrino mass in the 0.1 — 0.9 eV range. If this result is confirmed then it shows that 
neutrino masses are almost degenerate and well within reach of cosmological detection 
in the near future. 

Another important question which can be answered by cosmological observations is 
how large the total neutrino energy density is. Apart from the standard model prediction 
of three light neutrinos, such energy density can be either in the form of additional, 
sterile neutrino degrees of freedom, or a non-zero neutrino chemical potential. 

The paper is divided into sections in the following way: In section 2 I review the 
present cosmological data which can be used for analysis of neutrino physics. In section 
3 I discuss neutrino physics around the epoch of neutrino decoupling at a temperature of 
roughly 1 MeV, including the relation between neutrinos and Big Bang nucleosynthesis. 
Section 4 discusses neutrinos as dark matter particles, including mass constraints on 
light neutrinos, and sterile neutrino dark matter. Section 5 contains a relatively short 
review of neutrino physics in the very early universe from the perspective of leptogenesis. 
Finally, section 6 contains a discussion. 
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2. Cosmological data 

Large Scale Structure (LSS) - At present there are two large galaxy surveys of 
comparable size, the Sloan Digital Sky Survey (SDSS) [24, 25] and the 2dFGRS (2 degree 
Field Galaxy Redshift Survey) [26]. Once the SDSS is completed in 2005 it will be 
significantly larger and more accurate than the 2dFGRS. At present the two surveys 
are, however, comparable in precision. 

Cosmic Microwave Background (CMB) - The temperature fluctuations are conve- 
niently described in terms of the spherical harmonics power spectrum Ct,i = {\aim,\ 2 ), 
where ^(9,(p) = *}2i m ai m Yi m {6,(f)). Since Thomson scattering polarizes light, there are 
also power spectra coming from the polarization. The polarization can be divided into a 
curl-free (E) and a curl (B) component, yielding four independent power spectra: Ct,i, 
Ce,u C B ,u and the T-E cross-correlation C T e,i- 

The WMAP experiment has reported data only on Ct,i and Cte,i as described 
in Refs. [27-32]. Foreground contamination has already been subtracted from their 
published data. 

In addition to the WMAP experiment there are a number of other current CMB 
experiments, both ground and balloon based. Wang et al. [34] have provided a 
compilation of various data sets, and in addition to these there is the ACBAR experiment 
[33] which has measured the CMB at small scales. 

Other data - Apart from CMB and LSS data there are a number of other cosmological 
measurements of importance to neutrino cosmology. One is the measurement of the 
Hubble constant by the HST Hubble Key Project, H = 72 ± 8 km s" 1 Mpc" 1 [35]. 

The constraint on the matter density coming from measurements of distant type 
la supernovae is also important for neutrino physics. The most recent result is from 
the Supernova Cosmology Project [36] and yields Vt m = 0.25to'^l (statistical) ±0.04 
(identified systematics). 

3. Neutrino Decoupling 

3.1. Standard model 

In the standard model neutrinos interact via weak interactions with e + and e~. In the 
absence of oscillations neutrino decoupling can be followed via the Boltzmann equation 
for the single particle distribution function [37] 

~dt ~ ~dp = (1) 
where C co ii represents all elastic and inelastic interactions. In the standard model all 
these interactions are 2 <-> 2 interactions in which case the collision integral for process 
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i can be written 

1 f d 3 p 2 d 3 p 3 d 3 p 4 



Ccoll,i(/l) 



2E 1 J 2£ 2 (2vr) 3 2E 3 (2tt) 3 2E 4 (2tv) 3 
x (2n) 4 S\ Pl +p 2 -p 3 + P4 )A(f 1 J 2 J 3 J 4 )S\M\ 2 12 ^ 34 ., (2) 

where S\M\ 2 12 ^ 34i is the spin-summed and averaged matrix element including the 
symmetry factor S = 1/2 if there are identical particles in initial or final states. The 
phase-space factor is A(/i, / 2 , / 3 , / 4 ) = /s/^l - /i)(l - ^2) - /i/2(l - /a)(l - fa)- 

The matrix elements for all relevant processes can for instance be found in 
Ref. [38]. If Maxwell-Boltzmann statistics is used for all particles, and neutrinos are 
assumed to be in complete scattering equilbrium so that they can be represented by 
a single temperature, then the collision integral can be integrated to yield the average 
annihilation rate for a neutrino 

bGF (gl + g 2 R )T 5 , (3) 



7T 3 



where 
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This rate can then be compared with the Hubble expansion rate 

H=\.m g )J 2 ^ (5) 

to find the decoupling temperature from the criterion H = T\ t=Td . From this one 
finds that T D (v e ) ~ 2.4 MeV, T D (v^ T ) ~ 3.7 MeV, when ^ = 10.75, as is the case in 
the standard model. 

This means that neutrinos decouple at a temperature which is significantly higher 
than the electron mass. When e + e~ annihilation occurs around T ~ rn e /3, the 
neutrino temperature is unaffected whereas the photon temperature is heated by a 
factor (11/4) 1 / 3 . The relation T^/T 7 = (4/11) 1 / 3 ~ 0.71 holds to a precision of roughly 
one percent. The main correction comes from a slight heating of neutrinos by e + e _ 
annihilation, as well as finite temperature QED effects on the photon propagator [38- 
51]. 

3.2. Big Bang nucleosynthesis and the number of neutrino species 

Shortly after neutrino decoupling the weak interactions which keep neutrons and protons 
in statistical equilibrium freeze out. Again the criterion H = r| T=Tf can be applied 

to find that T frcczc ~ 0.5#* 1/6 MeV [37]. 

Eventually, at a temperature of roughly 0.2 MeV deuterium starts to form, and 
very quickly all free neutrons are processed into 4 He. The final helium abundance is 
therefore roughly given by 

2n n / rip 



1 + n n /n p 



(6) 

T~0.2 MeV 
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n n /n p is determined by its value at freeze out, roughly by the condition that 
n n /n p \ T =T tieme ~ e~( mn_mp )/ Tfrcczc . 

Since the freeze-out temperature is determined by g* this in turn means that g* 
can be inferred from a measurement of the helium abundance. However, since Yp is 
a function of both flbh 2 and g* it is necessary to use other measurements to constrain 
Q b h 2 in order to find a bound on g*. One customary method for doing this has been 
to use measurements of primordial deuterium to infer Q b h 2 and from that calculate a 
bound on Usually such bounds are expressed in terms of the equivalent number 
of neutrino species, N v = p/p uo , instead of g*. The exact value of the bound is quite 
uncertain because there are different and inconsistent measurements of the primordial 
helium abundance (see for instance Ref. [52] for a discussion of this issue). The most 
recent analyses are [52] where a value of 1.7 < N v < 3.0 (95% C.L.) was found and [53] 
which found the result N u = 2.5lo;g. The difference in these results can be attributed 
to different assumptions about uncertainties in the primordial helium abundance. 

Another interesting parameter which can be constrained by the same argument is 
the neutrino chemical potential, £„ = p u /T [54-57]. At first sight this looks like it is 
completely equivalent to constraining N u . However, this is not true because a chemical 
potential for electron neutrinos directly influences the n — p conversion rate. Therefore 
the bound on £ Ve from BBN alone is relatively stringent (—0.1 <^£j, e <^ 1 [54]) compared 
to that for muon and tau neutrinos ( <^7 [54]). However, as will be seen in the 
next section, neutrino oscillations have the effect of almost equilibrating the neutrino 
chemical potentials prior to BBN, completely changing this conclusion. 



3.3. The number of neutrino species - joint CMB and BBN analysis 

The BBN bound on the number of neutrino species presented in the previous section 
can be complemented by a similar bound from observations of the CMB and large scale 
structure. The CMB depends on N u mainly because of the early Integrated Sachs Wolfe 
effect which increases fluctuation power at scales slightly larger than the first acoustic 
peak. The large scale structure spectrum depends on N v because the scale of matter- 
radiation equality is changed by varying N u . 

Several recent papers have analyzed WMAP and 2dF data for bounds on N v 
[52, 53, 58-60], and some of the bounds are listed in Table 3.3. Recent analyses combining 
BBN, CMB, and large scale structure data can be found in [52,58], and these results 
are also listed in Table 3.3. 

Common for all the bounds is that N v = is ruled out by both BBN and CMB/LSS. 
This has the important consequence that the cosmological neutrino background has 
been positively detected, not only during the BBN epoch, but also much later, during 
structure formation. 
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Table 1. Various recent limits on the effective number of neutrino species, as well as 

t he data used. 

Ref. Bound on N v Data used 

Crotty et al. [59] 1.4 < N v < 6.8 CMB, LSS 

Hannestad [58] 0.9 < N v < 7.0 CMB, LSS 

Pierpaoli [60] 1.9 < N v < 6.62 CMB, LSS 

Barger et al. [52] 0.9 < N v < 8.3 CMB 

Hannestad [58] 2.3 < N v < 3.0 CMB, LSS, BBN 

Barger et al. [52] 1.7 < N v < 3.0 CMB, BBN 



3.4- The effect of oscillations 

In the previous section the one-particle distribution function, /, was used to describe 
neutrino evolution. However, for neutrinos the mass eigenstates are not equivalent to 
the flavour eigenstates because neutrinos are mixed. Therefore the evolution of the 
neutrino ensemble is not in general described by the three scalar functions, fi, but 
rather by the evolution of the neutrino density matrix, p = the diagonal elements 
of which correspond to f. 

For three-neutrino oscillations the formalism is quite complicated. However, the 
difference in Am 12 and Am 2 3, as well as the fact that sin2# 13 <C 1 means that the 
problem effectively reduces to a 2 x 2 oscillation problem in the standard model. A 
detailed account of the physics of neutrino oscillations in the early universe is outside 
the scope of the present paper, however an excellent and very thorough review can be 
found in Ref. [61] 

Without oscillations it is possible to compensate a very large chemical potential for 
muon and/or tau neutrinos with a small, negative electron neutrino chemical potential 
[54]. However, since neutrinos are almost maximally mixed a chemical potential in one 
flavour can be shared with other flavours, and the end result is that during BBN all 
three flavours have almost equal chemical potential. This in turn means that the bound 
on u e applies to all species so that [62-66]. 

161 = ^*0.15 (7) 

for i — e, /i, t. 

In models where sterile neutrinos are present even more remarkable oscillation 
phenomena can occur. However, I do not discuss this possibility further, except for 
the possibility of sterile neutrino warm dark matter, and instead refer to the review 
[61]. 

3.5. Low reheating temperature and neutrinos 

In most models of inflation the universe enters the normal, radiation dominated epoch 
at a reheating temperature, Trh, which is of order the electroweak scale or higher. 
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However, in principle it is possible that this reheating temperature is much lower, of 
order MeV. This possibility has been studied many times in the literature, and a very 
general bound of Trh ^ 1 MeV has been found [67-70] 

This very conservative bound comes from the fact that the light element abundances 
produced by big bang nucleosynthesis disagree with observations if the universe if matter 
dominated during BBN. However, a somewhat more stringent bound can be obtained by 
looking at neutrino thermalization during reheating. If a scalar particle is responsible for 
reheating then direct decay to neutrinos is suppressed because of the necessary helicity 
flip. This means that if the reheating temperature is too low neutrinos never thermalize. 
If this is the case then BBN predicts the wrong light element abundances. However, even 
if the heavy particle has a significant branching ratio into neutrinos there are problems 
with BBN. The reason is that neutrinos produced in decays are born with energies which 
are much higher than thermal. If the reheating temperature is too low then a population 
of high energy neutrinos will remain and also lead to conflict with observed light element 
abundances. A recent analysis showed that in general the reheating temperature cannot 
be below roughly 4 MeV [71]. 



4. Neutrino Dark Matter 



Neutrinos are a source of dark matter in the present day universe simply because they 
contribute to Q m . The present temperature of massless standard model neutrinos is 
T v Q = 1.95 K = 1.7 x 10~ 4 eV, and any neutrino with m ^> T v Q behaves like a standard 
non-relativistic dark matter particle. 

The present contribution to the matter density of N v neutrino species with standard 
weak interactions is given by 

^ = N w < 8 » 

Just from demanding that Q u < 1 one finds the bound [72, 73] 
46 eV 

4-1. The Tremaine-Gunn bound 

If neutrinos are the main source of dark matter, then they must also make up most of the 
galactic dark matter. However, neutrinos can only cluster in galaxies via energy loss due 
to gravitational relaxation since they do not suffer inelastic collisions. In distribution 
function language this corresponds to phase mixing of the distribution function [74]. By 
using the theorem that the phase-mixed or coarse grained distribution function must 
explicitly take values smaller than the maximum of the original distribution function 
one arrives at the condition 

/CG < /i/,max = ^ (10) 



Neutrinos in Cosmology 



8 



Because of this upper bound it is impossible to squeeze neutrino dark matter beyond a 
certain limit [74]. For the Milky Way this means that the neutrino mass must be larger 
than roughly 25 eV if neutrinos make up the dark matter. For irregular dwarf galaxies 
this limit increases to 100-300 eV [75,76], and means that standard model neutrinos 
cannot make up a dominant fraction of the dark matter. This bound is generally known 
as the Tremaine-Gunn bound. 

Note that this phase space argument is a purely classical argument, it is not related 
to the Pauli blocking principle for fermions (although, by using the Pauli principle 
f v <l one would arrive at a similar, but slightly weaker limit for neutrinos). In fact the 
Tremaine-Gunn bound works even for bosons if applied in a statistical sense [75], because 
even though there is no upper bound on the fine grained distribution function, only a 
very small number of particles reside at low momenta (unless there is a condensate). 
Therefore, although the exact value of the limit is model dependent, limit applies to any 
species that was once in thermal equilibrium. A notable counterexample is non-thermal 
axion dark matter which is produced directly into a condensate. 



4-2. Neutrino hot dark matter 



A much stronger upper bound on the neutrino mass than the one in Eq. (9) can be 
derived by noticing that the thermal history of neutrinos is very different from that of 
a WIMP because the neutrino only becomes non-relativistic very late. 

In an inhomogeneous universe the Boltzmann equation for a collisionless species is 

[77] 



L[f ] = RL = ^l + ^l^L + K^L = o 

Dt dr dr dx i dr dq l 



(11) 



where r is conformal time, dr = dt/a, and q l = ap % is comoving momentum. The second 
term on the right-hand side has to do with the velocity of the distribution in a given 
spatial point and the third term is the cosmo logical momentum redshift. 

Following Ma and Bertschinger [77] this can be rewritten as an equation for the 
perturbed part of / 



f{x\q\r)=h{q) [l + *{x\ q\ r) 
In synchronous gauge that equation is 



1 M dt! q xT dln/ 
jo or e dmq 



h + 6f] 9 

— 



= T C[f], 
Jo 



(12) 



(13) 



where qi 



qn J , fi 



n^kj, and e 



(q 2 



2U/2 



ti is the comoving wavevector. h 



and r] are the metric perturbations, defined from the perturbed space-time metric in 
synchronous gauge [77] 



ds 2 



hij 



[r)[-dT 2 + (Sij + hi^dx^da?], 



d 3 ke 



ih-x 



tikjhik, r) + (kikj - -5ij)6r)(k, r) 



(14) 
(15) 
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Expanding this in Legendre polynomials one arrives at a set of hierarchy equations 
4 2 • 

e =*■(!-,) 

'•' (iFn - (( + 1)F, +1 ) (16) 



2/ + 1 

For subhorizon scales (h — i) — 0) this reduces to the form 
S = - 4 




2<r = 

= 2TTI " (/ + (17) 

One should notice the similarity between this set of equations and the evolution 
hierarchy for spherical Bessel functions. Indeed the exact solution to the hierarchy is 

Fiikr) ~ jtikr) (18) 

This shows that the solution for 5 is an exponentially damped oscillation. On small 
scales, k > t, perturbations are erased. 

This in intuitively understandable in terms of free-streaming. Indeed the Bessel 
function solution comes from the fact that neutrinos are considered massless. In the limit 
of CDM the evolution hierarchy is truncated by the fact that = 0, so that the CDM 
perturbation equation is simply 5 = —h/2. For massless particles the free-streaming 
length is A = cr which is reflected in the solution to the Boltzmann hierarchy. Of course 
the solution only applies when neutrinos are strictly massless. Once T ~ m there is a 
smooth transition to the CDM solution. Therefore the final solution can be separated 
into two parts: 1) k > r(T = to): Neutrino perturbations are exponentially damped 
2) k < r(T = to): Neutrino perturbations follow the CDM perturbations. Calculating 
the free streaming wavenumber in a flat CDM cosmology leads to the simple numerical 
relation (applicable only for T cq 3> to 3> T ) [37] 



A 



FS 



20 Mpc ' ' " : / ° 7 '- 1 * 2xn 




(19) 



In Fig. 1 I have plotted transfer functions for various different neutrino masses in 
a flat ACDM universe (Q m + Vt v + Vt^ = 1). The parameters used were fi^ = 0.04, 
^cdm = 0.26 — VL V , Q a = 0.7, h = 0.7, and n — 1. 

When measuring fluctuations it is customary to use the power spectrum, P(k,r), 
defined as 

P(k,r) = \S\ 2 (r). (20) 
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0.01 0.10 
k (/i/Mpc) 



Figure 1. The transfer function T(k,t = to) for various different neutrino masses. 
The solid (black) line is for m u = 0, the long-dashed for m„ — 0.3 eV, and the dashed 
for m v = 1 eV. 



The power spectrum can be decomposed into a primordial part, Po(k), and a transfer 
function T(k, r), 

P(k,r) = P (k)T(k,r). (21) 

The transfer function at a particular time is found by solving the Boltzmann equation 
for S(t). 

At scales much smaller than the free-streaming scale the present matter power 
spectrum is suppressed roughly by the factor [78] 

AP(fc) = AT(fc,r = r ) ^ £2^ 

P(k) T(k,r = r ) { ] 

as long as Q u <C Q m . The numerical factor 8 is derived from a numerical solution of the 
Boltzmann equation, but the general structure of the equation is simple to understand. 
At scales smaller than the free-streaming scale the neutrino perturbations are washed 
out completely, leaving only perturbations in the non-relativistic matter (CDM and 
baryons). Therefore the relative suppression of power is proportional to the ratio of 
neutrino energy density to the overall matter density. Clearly the above relation only 
applies when Q u <C Q m , when fl„ becomes dominant the spectrum suppression becomes 
exponential as in the pure hot dark matter model. This effect is shown for different 
neutrino masses in Fig. 1. 

The effect of massive neutrinos on structure formation only applies to the scales 
below the free-streaming length. For neutrinos with masses of several eV the free- 
streaming scale is smaller than the scales which can be probed using present CMB data 
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Figure 2. Deviation of the best fit models for m u = (left) and 
Yj m v = 3 eV (right)- The CMB data shown are the WMAP binned data 
(http://lambda.gsfc.nasa.gov). The deviation is calculated as "Deviation in cr" = 

(-^model -^obs 



and therefore the power spectrum suppression can be seen only in large scale structure 
data. On the other hand, neutrinos of sub-eV mass behave almost like a relativistic 
neutrino species for CMB considerations. The main effect of a small neutrino mass on 
the CMB is that it leads to an enhanced early ISW effect. The reason is that the ratio of 
radiation to matter at recombination becomes larger because a sub-eV neutrino is still 
relativistic or semi-relativistic at recombination. With the WMAP data alone it is very 
difficult to constrain the neutrino mass, and to achieve a constraint which is competitive 
with current experimental bounds it is necessary to include LSS data from 2dF or SDSS. 
When this is done the bound becomes very strong, somewhere in the range of 1 eV for 
the sum of neutrino masses, depending on assumptions about priors. In Table 4.2 the 
present upper bound on the neutrino mass from various analyses is quoted, as well as 
the assumptions going into the derivation. 

As can be gauged from this table, a fairly robust bound on the sum of neutrino 
masses is at present somewhere around 1.0 eV, depending somewhat on specific priors 
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Table 2. Various recent limits on the neutrino mass from cosmology and the data 
sets used in deriving them. 1: WMAP data, 2: Other CMB data, 3: 2dF data, 4: 
Constraint on erg (different in 4 a and 4 b ), 5: SDSS data, 6: Constraint on Hq. 



Ref. 


Bound on m . 


„ Data used 


Spergel et al. (WMAP) [27] 


0.69 eV 


1,2,3,4«,6 


Hannestad [58] 


1.01 eV 


1,2,3,6 


Allen, Smith and Bridle [83] 


0.56i^ 6 eV 


l,2,3,4 b ,6 


Tegmark et al. (SDSS) [25] 


1.8 eV 


1,5 


Barger et al. [89] 


0.75 eV 


1,2,3,5,6 


Crotty, Lesgourgues and Pastor [81] 


1.0 (0.6) eV 


1,2,3,5 (6) 



and data sets used. 

It is also quite interesting to see what exactly provides this bound. It is often stated 
that the neutrino mass bound comes from large scale structure data, not from CMB 
because CMB probes larger scales. However, LSS data alone provides no limit on J2 m u 
because of degeneracies with other parameters (this is discussed in detail in Ref. [79]). 
On the other hand, WMAP in itself also does not provide a strong limit on the neutrino 
mass [25], because neutrino mass only has a limited effect on the scales probed by 
WMAP. Only the combination of the two types of data allows for a determination of 
J2 m u with any precision. Fig. 2 show deviation of the best fit models for J2 m u — 
and J2 m u — 3 eV from WMAP and SDSS data. From this figure it is obvious that 
models with high neutrino mass are not ruled out by any single data point, but rather 
by a general decrease in how well the combined data fits. One fairly evident problem 
with the high neutrino mass model is that the shape of the large scale structure power 
spectrum becomes wrong. The model spectrum has too much power at intermediate 
scales and too little at small scales. 

In the upper part of this figure the deviation of the best fit models from other 
cosmological data is shown. This data is not used in deriving the best fit models, and 
therefore the figure shows that the standard concordance model with J2 m u = is not 
only a better fit to CMB and LSS data, but also more consistent with other cosmological 
data. 

4-2.1. Combining measurements of m v and N v . The limits on neutrino masses 
discussed above apply only for neutrinos within the standard model, i.e. three light 
neutrinos with degenerate masses (if the sum is close to the upper bound). However, if 
there are additional neutrino species sharing the mass, or neutrinos have significant 
chemical potentials this bound is changed. Models with massive neutrinos have 
suppressed power at small scale, with suppression proportional to Q u /Q m . Adding 
relativistic energy further suppresses power at scales smaller than the horizon at matter- 
radiation equality. For the same matter density such a model would therefore be even 
more incompatible with data. However, if the matter density is increased together with 
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m U) and N v , excellent fits to the data can be obtained. This effect is shown in Fig. 3. 




Figure 3. Power spectra for ACDM models with fi& = 0.05, = 1, h = 0.7, n s = 1, 
and iV^massivo = 1 and a common large-scale normalization. The full line is for 0^ = 0, 
tt m = 0.25, N v = 3, dashed is for fi„ = 0.05, 9, m = 0.25, N v = 3, dotted is for 
VL V = 0.05, n m = 0.25, N v = 8, and long-dashed is for VL V = 0.05, VL m = 0.35, N u = 8 
(from [80]). 

The effect on likelihood contours for (Q, v , N v ) can be seen in Fig. 4 which is for the 
case where N v species the total mass equally. 

A thorough discussion of these models can be found in Refs. [80, 81]. 

4-2.2. Future neutrino mass measurements The present bound on the sum of neutrino 
masses is still much larger than the mass difference, Ar«23 ~ 0.05 eV [3,90], measured 
by atmospheric neutrino observatories and K2K . This means that if the sum of neutrino 
masses is anywhere close to saturating the bound then neutrino masses must the almost 
degenerate. The question is whether in the future it will be possible to measure masses 
which are of the order Am 23 , i.e. whether it can determined if neutrino masses are 
hierarchical. 

By combining future CMB data from the Planck satellite with a galaxy survey like 
the SDSS it has been estimated that neutrino masses as low as about 0.1 eV can be 
detected [84,85]. Another possibility is to use weak lensing of the CMB as a probe of 
neutrino mass. In this case it seems likely that a sensitivity below 0.1 eV can also be 
reached with CMB alone [92]. 

As noted in Ref. [85] the exact value of the sensitivity at this level depends both on 
whether the hierarchy is normal or inverted, and the exact value of the mass splittings. 
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Figure 4. Likelihood contours (68% and 95%) for the case of N v neutrinos with equal 
masses, calculated from WMAP and 2dF data (from [80]). 

4-3. Neutrino warm dark matter 

While CDM is defined as consisting of non-interacting particles which have essentially no 
free-streaming on any astronomically relevant scale, and HDM is defined by consisting of 
particles which become non-relativistic around matter radiation equality or later, warm 
dark matter is an intermediate. One of the simplest production mechanisms for warm 
dark matter is active-sterile neutrino oscillations in the early universe [104-108]. 

One possible benefit of warm dark matter is that it does have some free-streaming 
so that structure formation is suppressed on very small scales. This has been proposed 
as an explanation for the apparent discrepancy between observations of galaxies and 
numerical CDM structure formation simulations. In general simulations produce galaxy 
halos which have very steep inner density profiles p oc r a , where a ~ 1 — 1.5, and 
numerous subhalos [93, 94]. Neither of these characteristics are seen in observations and 
the explanation for this discrepancy remains an open question. If dark matter is warm 
instead of cold, with a free-streaming scale comparable to the size of a typical galaxy 
subhalo then the amount of substructure is suppressed, and possibly the central density 
profile is also flattened [95-101] . In both cases the mass of the dark matter particle 
should be around 1 keV [109, 110], assuming that it is thermally produced in the early 
universe. 

On the other hand, from measurements of the Lyman-o; forest flux power spectrum 
it has been possible to reconstruct the matter power spectrum on relatively small scales 
at high redshift. This spectrum does not show any evidence for suppression at sub- 
galaxy scales and has been used to put a lower bound on the mass of warm dark 
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matter particles of roughly 1.1 keV [102, 103]. An even more severe problem lies in the 
fact that star formation occurs relatively late in warm dark matter models because 
small scale structure is suppressed. This may be in conflict with the low-/ CMB 
temperature-polarization cross correlation measurement by WMAP which indicates very 
early reionization and therefore also early star formation. One recent investigation of 
this found warm dark matter to be inconsistent with WMAP for masses as high as 10 
keV [95]. 

The case for warm dark matter therefore seems quite marginal, although at present 
it is not definitively ruled out by any observations. 

4-4- Neutrinos as the source for high energy cosmic rays 

Nucleons with energy above the threshold for photo-pion production on the CMB rapidly 
downscatter in energy, the mean free path being of order 20-30 Mpc. At the present 
CMB temperature of 2.7 K the threshold, known as the Greisen-Zatzepin-Kuzmin energy 
[118], is roughly 4 x 10 19 GeV. On the other hand a significant number of particles 
with energies above the GZK energy have been observed. At present there is some 
controversy about the number of such particles observed between different experimental 
collaborations using different techniques. The HiRes [120] collaboration finds a decline 
in the number of events beyond the GZK energy which is in fact compatible with a 
cut-off. On the other hand the AGASA [119] collaboration finds that the spectrum 
is consistent with no cut-off, and even with a hardening of the spectrum at very high 
energies. Recent reviews of these issues can be found in Ref. [121]. 

Either these particles must come from relatively nearby sources or they are not 
nucleons (or nuclei) with standard model interactions. One explanation which has been 
proposed is the Z-burst scenario [111-115]. In this model, neutrino primaries with 
very high energy, E ^> Eqzk, annihilate with neutrinos in the cosmic background to 
produce the observed protons. The cross section is significantly enhanced at E CM = m z , 
corresponding to a neutrino primary energy of E ~ m 2 z /2E v $ ~ m 2 z /2m u . If neutrino 
masses are hierarchical then the largest mass is of order m v ~ Am atm ~ 0.05 eV, 
meaning that E v ~ 10 23 eV. If neutrinos are produced by pion decay in AGNs this must 
mean that particles of even higher energies are produced there. 

On the other hand, if neutrinos have masses close to saturating the cosmological 
bound then the primary energy can be significantly lower. 

Another requirement is that the neutrino annihilation must take place within the 
GZK sphere. This leads to a too low flux unless the rate is somehow enhanced. It has 
previously been proposed that this can be explained in models with significant neutrino 
chemical potential [116]. However, the present cosmological bounds on i] rule this out. 
If neutrinos are of eV mass then they do have significant clustering on GZK scales which 
can also enhance the rate by a factor of about 2. If ultra high energy cosmic rays are 
explained by the Z-burst this means that a mass bound on neutrinos can in principle 
be obtained. In Refs. [113, 114] it was estimated that if the annihilations are within the 
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galactic halo it requires a neutrino mass of m v = 2.34^1^1 eV, and if the annihilation 
happens within the local supercluster the mass must be m„ = 0.26±g;?2 eV. The first 
case is already ruled out, but the second possibility might work. It has, however, also 
been shown that the values obtained in [113, 114] are strongly model dependent [117]. 

At present the feasibility of the Z-burst scenario remains an open question. One 
problem is that the annihilation process produces a background of low energy gamma 
rays which may be in conflict with EGRET observations, depending on the magnitude 
of the local neutrino density enhancement. 

In any case the Z-burst scenario is also interesting from the possibility of getting 
an independent detection of the cosmological neutrino background if the needed very 
high fluxes of ultra-high energy neutrinos is measured in future detectors like Auger. 

5. Neutrinos in the very early universe - Leptogenesis 

A particularly attractive model for baryogenesis involves massive, right handed 
neutrinos, and is known as baryogenesis via leptogenesis [122]. The basic idea is that the 
masses of left handed Majorana neutrinos are generated from couplings to very massive, 
right handed neutrinos. 

These massive, right handed states are unstable and because they are Majorana 
particles their decay violate lepton number. Futhermore the decays are out of 
equilibrium and can violate CP which means that all the Sakharov conditions for 
generating a net lepton number are present. This lepton number can then subsequently 
be transferred to the baryon sector via standard model interactions and account for the 
observed baryon number of the universe. 

A particularly simple model for this is thermal leptogenesis where the right handed 
neutrinos are equilibrated at high temperatures directly via their interactions with the 
thermal plasma. In this case, the correct baryon number is produced only if the 
following conditions are fulfilled [123-128]: a) Masses of the light neutrinos must be 
less than about 0.1-0.15 eV, and b) Masses of the right handed neutrinos must be 
larger than about 10 8 GeV. The first condition is interesting because it provides a 
strong, albeit very model dependent, bound on the light neutrino masses. The second 
condition is interesting because it is so high that it might be in conflict with the upper 
bound on the reheating temperature in supersymmetric models. This bound arises from 
overproduction and subsequent decay of gravitinos and could probably be relaxed in 
models where the gravitino is the lightest supersymmetric particle. 

Taken at face value the thermal leptogenesis constraint on light neutrino masses is 
the most restrictive cosmological limit known. However, it is not a constraint at the 
same level as experimental bounds, or even bounds from CMB and large scale structure. 
The derivation involves a chain of assumptions: a) Leptogenesis is the correct model of 
baryogenesis, b) leptogenesis is thermal, c) The heavy neutrino masses are hierarchical. 
If either of the first two assumptions are relaxed then there is essentially no mass bound 
from this argument. If the last assumption is relaxed then it has been shown in Ref. [128] 
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Table 3. Summary table of cosmological neutrino mass limits. For completeness 
bounds from other sources, astrophysical and experimental, are also listed. 



Method 


Bound on J2 m v 


Data used 


0, /i 2 ^0.15 


14 eV 




CMB and LSS 


0.7-1 eV 


WMAP, 2dF, SDSS 


SN1987A 


m Ue < 5 - 20 eV (p e ) 


SN1987A cooling curve [129, 130] 


/3-decay 


m Ve < 2.2 eV (v e ) 


Mainz experiment [21] 


Oz/2/3-decay 


nT>v,ef£ < 0.35 eV 


Heidelberg-Moscow [131] 




0.1 eV < m^ eff < 0.9 eV 


Heidelberg- Moscow [22, 23] 



that the mass bound on the left handed neutrinos can be relaxed by almost an order of 
magnitude. 

With this in mind specific mass bounds on neutrino masses from thermal 
leptogenesis should be taken as both interesting and suggestive, but not as strict and 
generally applicable bounds. 

6. Discussion 

In the present paper I have discussed how cosmological observations can be used for 
probing fundamental properties of neutrinos which are not easily accessible in lab 
experiments. Particularly the measurement of absolute neutrino masses from CMB 
and large scale structure data has received significant attention over the past few years. 
In Table 3 I summarize neutrino mass bounds from cosmological observations and other 
astrophysical and experimental bounds. 

Another cornerstone of neutrino cosmology is the measurement of the total energy 
density in non-electromagnetically interacting particles. For many years Big Bang 
nucleosynthesis was the only probe of relativistic energy density, but with the advent of 
precision CMB and LSS data it has been possible to complement the BBN measurement. 
At present the cosmic neutrino background is seen in both BBN, CMB and LSS data 
at high significance. 

Finally, cosmology can also be used to probe the possibility of neutrino warm dark 
matter, which could be produced by active-sterile neutrino oscillations. 

In the coming years the steady stream of new observational data will continue, and 
the cosmological bounds on neutrino will improve accordingly. For instance, it has been 
estimated that with data from the upcoming Planck satellite it could be possible to 
measure neutrino masses as low as 0.1 eV. 

Certainly neutrino cosmology will continue to be a prospering field of research for 
the foreseeable future. 
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